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Summary. - The connection of unbroken SUSY quantum mechanics in its strictly 
isospectral form with the nonlinear Riccati superposition principle is pointed out. 
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m ' PACS 11.30.Pb - Super-symmetry. 

>' 

Generating quasi- isospectral potentials by the methods of supersymmetric quantum 
■ mechanics (SUSYQM) is a simple and well-known mathematical technique Q. Its 

relationships with intertwining, factorization, Darboux covariance, and inverse scat- 
tering procedures have been recently reviewed by the author g]. Also, the author, 
either in collaboration or alone || dealt extensively with a SUSYQM procedure 
based on the general Riccati solution that may be called the double Darboux 
general Riccati method (henceforth DDGR). For details the reader is directed to 
my previous works on DDGR and the references therein. The purpose of this note 
is to display the connection of the DDGR method with the Riccati nonlinear su- 
perposition principle. 

The DDGR procedure has been first used in physics by Mielnik for the quantum 
harmonic oscillator case ||. For simplicity reasons, I shall deal with the case of 
. unbroken SUSY, i.e., when one has at hand a true zero mode uo of the initial 

Schroedinger problem. DDGR leads to families of strictly isospectral potentials 
to the initial (bosonic) one. Asking for the most general superpotential (i.e., the 
general Riccati solution) such that the fermionic potential fulfills V+(x) = Wg + 
-gjr, it is obvious that one particular solution to this equation is w p , i.e., Witten's 
superpotential. Thus, one is led to consider the following Riccati equation w^ + 

= w p~^~~Sr' wnose general solution can be written down as w g (x) = w p (x)+^py, 
where v(x) is an unknown function. Using this ansatz, one obtains for the function 
v(x) the following Bernoulli equation 

^M-2v(x)w p (x) = l, (1) 



that has the solution 



To(x) + A 

v(x) = (2) 



where 2o( x ) — /* u o(y) dy, an d A is an integration constant considered as the free 
parameter of the method. Thus, w g (x) can be written as follows 



;(x; A) = w p (x) + In (2b (x) + A) = w p (x) + a (\) = - — 



where in (2) and (3) the well-known relationship (D = d/dx) 



In 



Uq(x) 
2 (x) + A 



w p (x) = -D In (u (x)) 



-<7 



(3) 



(4) 



has been used. The eg notation for the logarithmic derivative is borrowed from the 
book of Matveev and Salle |(| , and the subscript has been appended to show that 
the fixed solution on which the scheme is carried on is the Schroedinger ground 
state. Finally, one gets the family of "bosonic" potentials 



V-(x;A) 



2 dw g (x; A) 

W s (x;A) -^x— 

-v2r 



= V_(x)-2D 2 [ln(I (x) + A)] 
= V_(x)-2<7 , X (A) 

4u (xK(x) , 2 U 4(x) 



V_(x)- 



2 (x) + A (2 (x) + A) 2 



(5) 



All V_ (x; A) have the same supersymmetric partner potential V+ (x) obtained by 
deleting the ground state. By inspection of Eqs.(^) and (||) we can obtain the 
ground state wave functions for the potentials V_ (x; A) as follows 



u (x;A) = f(A) 



up(x) 
2o(x) + A : 



(6) 



where f(A) is a normalization factor of the form f(A) = yj A(A + 1). A connection 
with other isospectral methods has been found, by noticing that the limiting values 
-1 and for the parameter A lead to the Abraham-Moses procedure 0, and Pursey's 
one respectively. Here, we notice furthermore that by writing the parametric 
family in terms of their unique "fermionic" partner 



V_(x;A) = V+(x)-2D 2 



In 



2 (x) + A 



no 



V+(x) - 2D 2 In 



u (x; A) / ' 



(7) 



the inverse character of the parametric Darboux transformation is manifest Q. 
Therefore, a two step interpretation is possible. In the first step, one goes to the 
fermionic partner system, and in the second one returns to a A-distorted bosonic 
system. This is why it can be called a double Darboux transformation, although 
this may be somewhat confusing because of Adler's method ||, which is also a two 
step construction and is known by the same name. The strictly isospectral DDGR 
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method can be applied to any one-dimensional system whose dynamics is dictated 
by a Schroedinger-like equation. Moreover, one can employ combinations of any 
pairs of Abraham-Moses procedure, Pursey's one, and the Darboux one. However, 
only the DDGR leads to reflection and transmission amplitudes identical to those 
of the original potential. 

We now make use of the fact that Riccati equations are differential equations 
admitting a nonlinear superposition principle |l0[ , i.e., it is possible to write down 
the general solution as a combination of particular solutions without using any 
quadrature. This is due to the well-known property JHJ that given three particular 
solutions one can obtain the general Riccati solution from the following k invariant 

w - Wl w 3 - Wl 

: = k, (8 

w — w 2 w 3 — w 2 

as the superposition 

kwi(w 3 - w 2 ) + w 2 (wi - w 3 ) 

w — T7 r- • l y j 

k(w 3 — w 2 j + wi — w 3 

In particular, we have seen that within DDGR the general Riccati solution can be 
written in the form w(A) = — do + uo(X). Considering three particular solutions 
w(Aj), i=l,2,3, one can write the above invariant as 

<7o(A) - g (Ai) _ g (A 3 ) - g (Ai) = k ^ iQ ^ 
cro(A) - ct (A 2 ) cr (A 3 ) - cto(A 2 ) 

i.e., in terms of the A-dependent logarithmic derivatives. The nonlinear Riccati 
invariant k can be seen as a constraint on the Darboux contributions to the para- 
metric bosonic potentials or as a constraint on a sequence of their free parameters 
Ai if written as follows 

/ Qq,x(A) - J <tq,x(Ai) _ / o- ,x(A 3 ) - / <tq,x(Ai) = k 
/ Co,x(A) - / cr ,x(A 2 ) ' / cr 0iX (A 3 ) - / cr , x (A 2 ) 

It is worthwhile to notice that an even more general form of the Riccati k 
invariant can be written down if one do not stop at the one-parameter type of 
bosonic Riccati solution. In another work Jl2[ , I have provided a simple multiple- 
parameter (iterative) generalization of Mielnik's construction. At an arbitrary i- 
order of iteration, one can write the general Riccati solution as Wg 1 ' = w + w\ ± + 
...~w\ i l -fwA, where up to i— 1 the A parameters are fixed and only the last parameter 
is a free one, whereas wv = Dln(Aj + f x Fj_i). The F's are integration factors of 

the type Fj = e~ I 2w p' i where Wp' is the particular Riccati solution at order j, 
which is of the same form as Wg' but with no free parameter. Since all these Riccati 
solutions correspond to the same fermionic Schroedinger potential, any triplet of 
multiple-parameter particular solutions that may be of arbitrary hierarchical order 
j can enter in the Riccati k invariant and may be used to get the general solution 
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at an arbitrary order. Thus, one can write 

/a»(..,A)-/^(...,A (jl) ) /^ ) (-,A j3 )-/<T ^(...,A jl ) k 

Finally, let us notice that instead of the zero mode solution one can employ 
any other solution of the initial Hamiltonian Hbu(x) = eu(x), where the so-called 
factorization energy e should be less than the ground state energy of Hb, and an 
analytic continuation is performed in the independent variable |i~3| to turn u into a 
nodeless function. In this way no singularities with respect to the initial potential 
are introduced. For this more general case the zero subindex in all the above 
equations is to be droped. 
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